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We employ the efTective field theory approach for multi-field inflation which is a generalization 
of Weinberg's work 1]. In this method the first correction terms in addition to standard terms in 
the Lagrangian have been considered. These terms contain up to the fourth derivative of the fields 
including the scalar field and the metric. The results show the possible shapes of the interaction 
terms resulting eventually in non-Gaussianity in a general formalism. In addition generally the speed 
of sound is different but almost unity. Since in this method the adiabatic mode is not discriminated 
initially so we define the adiabatic as well as entropy modes for a specific two-field model. It has been 
shown that the non-Gaussianity of the adiabatic mode and the entropy mode are correlated in shape 
and amplitude. It is shown that even for speed close to unity large non-Gaussianities are possible in 
multi-field case. The amount of the non-Gaussianity depends on the curvature of the classical path 
in the phase-space in the Hubble unit such that it is large for the large curvature. In addition it is 
emphasized that the time derivative of adiabatic and entropy perturbations do not transform due 
Cn ' to the shift symmetry as well as the original perturbations. Though two specific combinations of 

them are invariant under such a symmetry and these combinations should be employed to construct 
an effective field theory of multi-field inflation. 



o 

(N 



^H I PACS numbers: 98.80.Cq 



Contents 
I. Introduction 



fS |. II. Briefly Review of Weinberg's Approach [l|] y 
D 

^ III. EfTective Field Theory of Multi-Field Inflation [3 

^__l IV. A Speciflc Case: Adiabatic versus Entropy Perturbation 5 

J> , A. The Most General Two-Field Model 5 

\^ ■ B. Adiabatic vs. Entropy Modes 6 

^O ' 1. The Amplitude of Non-Gaussianity 8 

CN ! 2. The Shape of Non-Gaussianity IC 

"^ _ C. An Example IC 

CO . D. Some Clarifications on Differences with Senatore and Zaldarriaga [J] [ll 

O ^ 

CN ■ V. Conclusions m 



Acknowledgments [ij 

/\^ , A. The Most General Lagrangian for Multi-Field Inflation \ll 
H ■ 

.'?.■ B. Comparison With Gordon et al. (Tol| ^ 

References Ul 



'Electronic address: l mma@aims.ac.zal 



I. INTRODUCTION 

It is generally believed that inflation can be a solution to the problems of standard cosmology such as the horizon, 
flatness and monopole problem. In addition to these achievements, inflation's predictions are compatible with the 
large scale structure and CMB fluctuations which is strong evidence in favour of inflation. The idea of inflation is the 
existence of an exponentially expanding universe at early times. But identifying a unique theoretical realization of 
this period is challenging. Many theoretical models are compatible with the observational data. For example, they 
are in agreement with adiabatic, nearly Gaussian fluctuations in the CMB fluctuations. To potentially discriminate 
between them more accurate observations, such as PLANCK, are needed. This fact is a starting point for a huge 
amount of work on studying non-Gaussianity of primordial fluctuations. In this field the effective field theory approach 
to inflation has been used to study the general possible interaction terms in the single field models [i|-|3| and in the 
multi-field context [J, Q. 

The advantage of using effective theories can be seen in two regimes. Sometimes a full theory exists for an energy 
domain of interest. In this case the effective field theory may be performed to simplify calculations in a special sub- 
domain of energy. In the second case the full theory is not known for the energy scales of interest. Here, by imposing 
the symmetries of the full theory one can still build an effective field theory. In this situation the most general form of 
the allowed theory, e.g. a general Lagrangian, is constructed; by comparison to observations unspecified coefficients 
can be fixed. Eventually the deduced effective theory may shed some lights on the real theory which is beyond our 
current understanding. In the special case of infiation in addition to the above reasons, the effective field theory 
approach can be used to justify the use of scalar fields as inflatons, as well as to provide a systematic classification of 
non-Gaussianities [6| among other properties. 

In ^3| the effective field theory has been developed for an inflationary single field model. In their approach the 
Lagrangian is determined by all spatially diffeomorphism invariant operators. Then the broken time invariance is 
reproduced by a scalar field which transforms in a definite form under diffeomorphism transformation. This scalar 
field is well-known as the Goldstone boson. It is shown that this scalar field represents the curvature perturbation 
in the validity regime of the effective field theory. In [J] the generalization to multi-field inflation is studied. The 
existence of more than one field in the early universe is not unnatural and the extra fields may have observable 
consequences. For example entropy modes (a property of multi-field models) can affect the curvature mode which is 
for example in the CMB. In [l| an alternative approach to Q has been given for the effective field theory of single 
field inflation. In this approach all the possible terms containing up to the fourth derivative of a scalar field and the 
metric enter the effective Lagrangian. The final result is in agreement with , except some additional fourth ordered 
contributions tracing back to geometrical terms. In this approach, due to the presence of metric perturbations, it is 
possible to study the gravitational wave behavior which differs in the propagation of waves with different helicities. 

In this work we are going to generalize Weinberg's approach |lj to multi-field models. In the following we will avoid 
the scalar metric perturbations by an appropriate gauge choice. However it is mentioned that for energy scales of 
interest the existence of them has no observable effects on non-Gaussianity [a, 0I- Note that in [ll, [2:l3i the additional 
correction terms arise via space-time derivatives of the perturbations. However one can extend the effective field theory 
formalism to include correction terms corresponding to the potential terms. It is mentioned in [8| that they have no 
significant contribution to non-Gaussianities since they are highly restricted by the effectiveness of the inflationary 
era. But it is well-known that in the context of multi-field inflation the non-Gaussianity window becomes wider and 
maybe observable by the future data ^Qj . This fact also has been considered in the context of effective field theory for 
multi-field inflation in different aspects Q. 

In the next section we briefly review the main results of p|. In the third section we generalize the idea of [l] to 
illustrate the perturbations in the most general multi-fleld model. This section is based on the first appendix where 
we find the most Lagrangian for multi-field models. Then in the fourth section we concentrate on a two-field case, 
studying the evolution of adiabatic and entropy modes in details. In this section we will discuss on the amplitude and 
shape of non-Gaussianity in our model and illustrate a specific example. At the end of this section we infer to some 
differences between this approach and Senatore and Zaldarriaga [4]. In the second appendix we compare our results 
with jlfll as a check. Finally we conclude in the last section. 

II. BRIEFLY REVIEW OF WEINBERG'S APPROACH [ij] 

To generalize Einstein-Hilbert action in the presence of matter field it is possible to add the terms containing higher 
order derivatives in the Lagrangian in addition to the standard second order ones. In principle these additional terms 
can be important in the larger energy scales. As discussed in ^1] this situation occurs naturally in the inflation era 
before horizon crossing. According to the observations; the Hubble parameter, iJ, and physical momentum, fc/a, are 
equal (at horizon crossing) and much less than Mp and even unification scale. But due to denominator of physical 



momentum in a period before horizon crossing the physical momentum has had larger value. As a consequence, 
considering the correction terms will help us to understand better the inflationary predictions. 

Due to the above discussions Weinberg in [l| has studied the effects of the fourth order derivative terms in the 
Einstein-Hilbert Lagrangian in the presence of one scalar field. We are going to generalize this model by adding more 
than one scalar field which is interesting for its well-known observational consequences. Before that let us review very 
briefly [i^l the main results of [l]. The starting point is the Einstein-Hilbert Lagrangian which includes the leading 
term 
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where dimensionless ip = tpc/M is defined such that the kinetic term of ifc has the canonical form. Obviously Lpc has 
dimension of mass. It is now easier to define the hierarchy of different derivative terms as an advantage of introduction 
of the scale M explicitly [ig. The leading correction terms are satisfied general covariance and contain four derivatives. 
These term can be reduced to the following form 



^C = V5/(^) ("g^V,^^,.") + ^h,{ip)C^^-'"'C^,p 



^gh2{v)e'^^'^'^Cj''C,,p„ 
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where /, hi and /12 are some dimensionless arbitrary functions which are assumed to be in order onefly] and C^upa 
is the Weyl tensor. It is noteworthy that the above terms are not all the generally covariant terms containing four 
derivatives. But all the allowed terms except the above terms are transformed to these terms by employing the 
equation of motion for the leading term as well as ignoring the surface terms (for details see Appendix |X]). For the 
scalar perturbations it is convenient to assume a gauge in which metric scalar perturbations vanish. In this gauge by 
splitting the scalar field to its background and perturbed parts as Lp = (p + 5ip the Lagrangian becomes 



-C = V5 
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+ 4a3/((^)<^2 L-^^ ^ ^-^3/^ _ a'^5^{V5^f/^ + ^^'//^' - ]-a-^5^^ {V5^f /^^ + '^a'^i^S^f /^^ 



which reduces to the Lagrangian (19) in [l[ with tt = Sip/^ and H — —Cp^{]\''P -I- Af{(p)^'^)/2Mp. This result 
is compatible with [2] with a minor disagreement. This disagreement is in the presence of quartic terms as well as 
quadratic and cubic terms. Due to the above Lagrangian obviously ignoring the correction term, i.e. setting /((/?) = 0, 
results in a model with c^ = 1, where Cs is the speed of sound. But in the presence of the correction term the speed 
of sound is not one and may cause large non-Gaussianity. In addition the terms in the third line of ^ infer to the 
possible shapes of non-Gaussianities as well as their amplitude. 

In this section we very briefly reviewed the idea of [l| for scalar perturbations in the context of effective field theory 
for inflation. In addition to scalar perturbation in jl| the tensor perturbations have been considered. It is concluded in 
[ij] that the propagation of gravitational wave depends on the helicity of the wave in this model. In the next section we 
are going to generalize the above idea for a multi-fleld theory of inflation without considering the tensor perturbations. 
Since existence of multi-scalar-fleld has no effect on the tensor perturbations and consequently gravitational wave. 
The detailed calculations for the next section is in appendix |X] which is also useful for clarifying the case of one field 
studying very briefly in this section. 



III. EFFECTIVE FIELD THEORY OF MULTI-FIELD INFLATION 

The corresponding Lagrangian to ^ for multi-fleld inflation can be written as follow, which has been deduced in 
details in the appendix 1X1 
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which exactly reduces to ^ for one field case[20|- Now, by splitting the scalar fields to their background and perturbed 
parts ifi = ifj + Sip we are going to study the Lagrangian for the perturbations as well as the background. To do 
this we start with ^ without worrying about the tensor perturbations. The above Lagrangian can be written as 
C = Cf) + A£ such that 
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where a = a{t) is the scale factor of the FRW metric and di are spatial derivatives. Note that the terms containing 5ipi 
without any differentiations do not show themselves in the Lagrangian effectively. Since the n*^ equation of motion 
causes vanishing of the coefficients of (n + 1)*'* terms without any differentiation. Also the linear perturbation terms 
even including differentiation vanish because of the same reason. So effectively the Lagrangian for the perturbations 
is 
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The second, third and fourth order of perturbations respectively can be represented as follows 

a"^A/:(^) = hl-'^^{if>)l 0j0jS(pi^6ipL + 0i0K5(pjSip^ + 0j0LSipj^Sipj + 0K0jSfjSipi^ + 0L0,jSipjSip]^ + 0K0LSipjSLpj 
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a ^A£^^^ = b"^^{ip)<0iS(pjSipj^6ifj^ + 0jS^pjSipj^Sipi^ + 0KS(pLSipjSipj + ipLSipjiSipjSipj (6) 

- a~^ ( 0iSipjdiSipKd''5ipL + 0jSipjdiSLpKd'-SipL + 0KSipLdi5ipid'Stpj + 0LSipxdiSipid'-5ipj 
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Note that the above result exactly reduces to single field result in ^ with hl^'^^^{(p) = /((^)[2lj. It is obvious from 
dS]) that the speed of sound is not one in the presence of non- vanishing bl'^^'^Up). Note that due to b^''^^{(p) the 
cubic and quartic terms are appeared. This feature is in disagreement with [2, lj|- In their work the coefficient which 
displays Cs just connects to the cubic term. But here it connects to the fourth order term too. In the next section we 
restrict the model to a two-field model. This makes it possible to study the adiabatic and entropy perturbations in 
more details without loss of generality in the main results. 

IV. A SPECIFIC CASE: ADIABATIC VERSUS ENTROPY PERTURBATION 

In this section we re-write the above formalism in the language of adiabatic and entropy perturbations for a two- 
field model. This is crucial in this approach since in contrast to Q here the adiabatic perturbation is not initially 
supposed. In [J| the additional perturbations are added to a model already containing the adiabatic perturbation 
i.e. [2|. In Q the Goldstone boson, introduced in [2|, plays the role of the adiabatic perturbation and the additional 
fields are employed as the entropy perturbations. But in our model there is no initially difference between (p/'s and 
consequently Sfi^s. So it is critical to distinguish between adiabatic and entropy modes to manifest their different 
physical meanings. 

A. The Most General Two-Field Model 

In this subsection we re-do perturbation calculations for a two-field model. To do this we start with (j3]) for two 
fields named (p and x 

/: = -«'! - ^d,ipd^ip-^d^xd^x-Ml.U{ip,x)+9i{^,x){d,ipd^^f+g2iip,x){d^xd'xf (8) 
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+ 95{^,x){d^vd''^){d,xd''x) + 96{^,x)id^vd''x){d.^d''x) 

where a — a{t) is the scale factor and g^'s are some arbitrary dimensionless and order one functions as mentioned 
before. By assuming ip ^ (f + Sip and X — X + ^X the above Lagrangian reduces to 

a-'Co = (-^d^^d'^^-^d,xd''x-Ml,U{^,x)+gi{d^^d^^f+92{d^xd^x)' (9) 
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= —'f'^ + —X^ ~ MpU{<fi, x) + 9iV^ + 92X^ + 53<^^X + 9WX^ + (.95 + 9&Wx^ 
for the background part. The corresponding equations of motion for (p reads as 
d 
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where ' is the differentiation with respect to (p and the similar equation is true for x- It is straightforward but messy 
to show that the above equation of motion (as well as x's) for the background causes the Lagrangian of the first order 
perturbation becomes vanishing. So the non-trivial terms start from the second order perturbations succeeding with 
the third and the fourth order terms [22l for ([5]) as the following 
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It is interesting to mention that for ip — x^ i-e. going back to one field case, all the above relations reduce to ^ with 
f = 91 + 92 + 93 + 34 + 95 + .96 as it was expected. 
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B. Adiabatic vs. Entropy Modes 

Now let us re-write the above terms in the language of adiabatic and entropy perturbations. The adiabatic pertur- 
bation is along the classical path and the entropy perturbation is orthogonal to it. Due to jlO| they can be defined as 
follows 

Sa = f.S, Ss = N.S (13) 

where Sa and Ss are the adiabatic and entropy modes respectively and 

S=iSip,6x), T = (cos 6*, sin 6^) = (ip/a, xja) , 

where a^ = ip'^ + x^. One can show easily 

Sip = cos 9 f .6 + sine N.S, 
and it is easy to see that 

Sa = eSs + f.S, 

and due to above relations 

diSip = cos9 diSa + sinO diSs, diSx = sm6 diS a — cos 6 diSs. 

Now by the above definitions we re-write the results of previous sub-section by plugging Sa and Ss in. Let's start 
with the kinetic terms for leading order term in (flO)) 
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the same procedure is applicable for a ^ ( -^diSipd'^Sip H — 2^9i(5x(9Mx) in (ITUl) . Now let us assume Afi = M2 — M 
to make it comparable to results in [10]. In this case 
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According to the above Lagrangian the equations of motion for Sa and 6s become [23| 

da + 3HSa - a~^d'd^5a + (K<t - d^)Sa ^ 2eSs + (6 + iHO - V„s)5s 
5s + iHSs - a-'^d''d^5s + {Vss - 0'^)5s = -29da - [9 + 3iJ^ + V„s)5a 



(17) 



where Kts = (cos^ 6* — s\v? 9)V^^ + sin0cos0(Fx;^ — ^^'v)- The above results are exactly same as (47) and (48) in 
[lOJ I when ignoring metric perturbations, see Appendix B. Now let us do the same procedure for the second order 
perturbations due to the first order correction term in (jlOp . At the first the terms containing time derivative 
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which can be re-written as the following j24 
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gi cos^ 9 + g2 sin^ 9 + g^ cos^ 6* sin + 54 cos 6* sin"^ ^ + (55 + ge) cos^ 9 sin^ 
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3a^{f.t){N.^ X 
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and similarly for the spatial differentiation 
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Up to now we fully considered the second order perturbation terms in the language of adiabatic and entropy modes. 
In quadratic level the speed of sound is a matter of interest hence it is noteworthy to take a look at it. It is obvious to 
the above relations that Sa and Ss can have different speeds of sound generally. To have a sense about it let's consider 
a special case that 6a has Cs — I and Ss has Cg ^ 1. To see this, assume the special case with gi = g2 = 93 = gi = 
and 35 + 36 = 0. In this case the coefficients of Sa and (diSa)'^ are same and result in Cg = 1 for 5cr [25|. But the 

coefficient of Ss is ^ and for (diSs)^ is —^a~'^{M^ — 2g<^cT^) that means c^ — 1 — 2^j^. Note that here we write the 
M explicitly to make the comparison of the terms easier. A characteristic property of c^ is its a^-dependence which 
seems interesting. However the effective field theory is valid where the correction terms are smaller than the leading 

\I,IJKL /-^,-2 IJ,-f-f^i/ 7;M^.2 
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< 1. Even more, 
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<< 1 should be satisfied 



terms in Q to have an acceptable expansion i.e. 

to make skipping higher order correction terms in ^ acceptable. So the speed of sound in this model is almost one. 
This fact shows that for the single field model the large non-Gaussinity is not expected. But in the following we will 
discuss on the case of multi-field models. In multi-field models even with c^ ~ 1 the large non-Gaussianity can be 
occurred in some specific circumstances. 

One can do this procedure for the higher order perturbation terms ([TlT) and (TT^ which are the fundaments of 
studying non-Gaussianity for multi-field models. The results are as follows[26] for a^'^l\0-'^> 
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+ ct(7V.J)^ X 

- 34 COS* 9 + g^ sin"* — 21 2^2 ^ (55 + ffe) ) cos'^ sin 6* -I- 2 1 2gi — (35 + ^e) ) cos 9 svc? 9 + 3{g4 — gs) cos^ 9 sin^ 9 
3CT(f.i)^(7V.i) X 

- 33 cos" 9 + g4: sin"* 9 + 2(2gi - {g^ + ge)] cos^ 6* sin 6* - 2 ( 2^2 - (35 -H ^e) ) cos 6* sin^ 6* -t- 3(33 - 34) cos^ 9 sin^ 9 
2&{f.^{N.t)^ X 

(.95 + .96) ( cos"* 9 + sin"* 9] + 8(34 - 33) ( cos^ 6* sin 6* - cos 9 sin^ ^ ) + ( Hqi + .92) - 2(35 + 36) ) cos^ 9 sin^ 



and for a ■^AC^'^'> it becomes 

{T.6) X gi CDS'* 6 + g2 sin^ 9 + g^ cos^ 9 sin 6 + g^ cos 9 sin^ 6* + (55 + gg) cos^ 9 sin^ 6* 

+ (N-S) X 52 cos"* + gi sin'' — 94 008'^ 0sinf? — 53 cos6'sin^ 6* + (95 + gg) cos^ 0sin^ 

+ (f.(i)^(7v.i) X 

- 33 CDS'* 9 + g4 sin'' 6* + 2 f 2^1 - (55 + ge) ) cos^ 61 sin 6* - 2 f 2^2 - (ffs + 55) ) cos 61 sin^ 9 + 3{g3 - g^) cos^ 61 sin^ 9 
+ (f.l)(iV.i)^x 

- 34 cos" 9 + g3 sin" 6* - 2 ( 2.92 - (.95 + ffe) ) cos^ 6* sin 6* + 2 ( 2gi - (^5 + c 



4o- 



cos 9 sin^ ^ + 3(54 - 53) cos^ 6* sin^ 9 

{95 + ge) ( cos" 9 + sin" 9 J + 8(^4 - gs) ( cos^ 6* sin 6* - cos 9 sin^ ^ ) + 2 ( 3(5i + 52) - 2(35 + gg) ) cos^ 9 sin^ 6* 
Now let us just consider the terms containing (T.S) and write them together as[27f 
gi cos 9 + g2 sin 9 + 93 cos^ 9 sin 9 + g^ cos sin'^ + (95 + 90) cos^ 6* sin^ 



x((f.^y + (f.jy/a + i(f.5)V 



Comparison the above relation with the relation in (|3]) manifests that a, [...] and (T.S) play the role of (p, f{(p) and (5(p 

respectively. The significant property of this model is the appearance of {T.6j and {N.dj or equivalently {da — 96s) 
and (5s + 9da) respectively. This means that Sa and 6s are always together and the same for Ss and Sa. This 
characteristic feature of this model has some observational consequences which are discussed in the following. 



1. The Amplitude of Non-Gaussianity 

Now we are going to estimate the non-Gaussianity amplitude. To do this one procedure is comparison between the 
non- linear terms and the linear ones. Mathematically, the amplitude of non-Gaussianity /jvl, bi-spectrum, can be 
estimated as ^j^ x C~^ where C, is the curvature perturbation Q. Note that the dominant amplitude of the terms 
containing time derivatives comes from their amplitude at horizon crossing. At this time 4; ^ H where H is the 
Hubble constant. Hence for the second order perturbations the Lagrangian C^'^' can be written in an abstract form 
as {H^,H9,9^} x A'P x Sa^. The same analyze for £(3) results in {H^,H^9,H9^,9^} x f{g^) x ax Sa^. So in an 
abstract form 



£(3) {H^,H^9,H9^J^} ( f{g, 



^^!^2 ) X ^ ^ {H\H^9,H9\9^) ^ ( fj^^^ , ^ ^ 



(20) 



/:(2) {H^,H9,9^} \ M^ J & Hx{m,H9,9^} \ M^ 

where C ~ -^^ is interpreted as curvature perturbation. Now we consider two different regimes 9 « H and 9 » H. 
The first regime, 9 << H, physically means that the model is a single field model effectively. In this case the amplitude 
of bi-spectrum can be approximated by 



Tnl 



/:(^V_i fig.)^: 



£(2) 



IVP 



(21) 



but remember that the validity of the effective field theory imposes JI2 o"^ < 1- So in this case as mentioned before 
there is no significant non-Gaussianity which is in agreement with the single field models of inflation [11| . But the 
other case, 9 » H, means that the classical path in the phase space is highly curved [l^l- In other words it means 
the classical path in the phase-space is far from a straight line {9 = 0). So the existence of the entropic held is 
unavoidable. For this case the amplitude of /atl is 






(22) 
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Now the factor jj x ( |f:l o-^) can be large and results in large non-Gaussianity consequently. So the large curvature 
of the classical path in the phase-space results in the large non-Gaussianity. Though this result can be compared to 
the other works in the literature [3] and in the effective field theory context Q but the curvature of the classical path 
is restricted due to observed scale invariant power spectrum [l3j . 

For a moment let us relax the constraint on the j^i (y^ ■ Consequently the correction terms in ([5|) causes the large 
non-Gaussinity. However the relaxation of the constraint can be justified by assuming that our model is completely 
described by (|S]) without any higher order correction terms. This needs fine tuning which is not impossible but it 
is not natural. However there is another method to rationalize this assumption. Instead of fine tuning the model 
automatically shows this property via for example Vainshtein mechanism |14l . |l5| . 

The same is applicable for tri-spectrum by an estimation as t^l ^ Tijj ^ C~^- The fourth order Lagrangian, C'^^\ 
can be written in the abstract form as {H'^,H^d,H'^e^, 116^,119^,6'^} x f{g^) x Sa^ and then 



TNL ^ T^T^C ~ ^rrrf^ (23) 



ioi 6 « H and 



for 6 » H. 



^^'K-2 / M .. ffi9^>■2 



--Z(^C-H^ X ^.^). (24) 



2. The Shape of Non- Gaussiamty 

Now let us focus on the shape of possible non-Gaussinity predicted by our model. In principle all the possible 
interaction terms between 5a, 5s and their derivatives exist in our model. This fact means all the non-Gaussianity 
shape can be produced. However the amplitude of different shapes are different. As a general argument it can be 
emphasized that for different regimes of -^ different shapes are dominant. For the case 6 « H the terms containing 
time derivatives become dominated. This means in this limit the equilateral shape is the main one among the others. 
But it does not mean the other shapes do not exist i.e. the "Cosine" between different shapes are not zero. In the 
other limit, 6 » H, the terms without derivative become dominant and then the local shape is dominant. This result 
is in agreement with the prediction for multi-field inflation models [llj. Note that in [J] since the additional entropy 
perturbations are added by symmetry then they do not have any term without derivative in their Lagrangian. In this 
sense they do not predict a dominant local shape for their model which is in disagreement with our result. 

A characteristic feature of this model is appearance of just two combinations of the fields i.e. T.5 and N.6 in all 
the terms including the second, third and fourth orders. To explain what is the physical result of this fact let us 

concentrate on T.5 = Sa — dSs, as an example. The third order term of this combination is {T.5)^ = Sa^ — W5a'^5s + 
W^5a5s^ — 9^5s^ . Without worrying about the amplitude in this part let us focus on the first and the last term. The 
definite prediction of this model is that if any equilateral non-Gaussianity due to the first term, i.e. 5(7^ , be observed 
then it has to be observed a local non-Gaussianity due to the last term[23|, i.e. 5s^ . So the non-Gaussianity predicted 
by this model cannot be pure e.g. pure equilateral shape. Hence mathematically, the "Cosine" between two shapes 
cannot be zero and more the "Cosine" depends on the 6 and is fixed by the model. This argument is true for the other 
third order terms as well as fourth order ones. To conclude, it seems that this model predicts a definite combination 
of different shapes for the non-Gaussianity if the amplitude allows to observe them. 

C. An Example 

In this subsection we are going to show how the general statements mentioned before do work in a simple example. 
Here we assume that all the giS vanish except gi{ip, x)12MI which is a generalization of the model in [3]. In addition 
we assume there is no potential term[30|. According to the background Lagrangian ^ the equations of motion for 
our special case become 

31 -2^ , orr -/l , . 31 -2^ , 1 ■ 3 /^o • ^.9l , a ■ ^9l 
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where Af — Mi = M2 is assumed. On the other hand, what can cause the significant non-Gaussianity is 6 as 
mentioned before. In general due to the definition of 6 in ([T?l) . 6 can be read as 



(p2. _)_ ^2 



and in our special case by considering (j25p it becomes (up to the first order of gi -^ 



-^24(-5i-^'^ ^^ ' ' ' ^' 



H \ "'iVPjip^+x^ HA'Pip^+x^ 



C'(y9 OX 



2^ 



The condition ffi-p^ < 1 ensures the validity of the effective field theory. So it is not bad to estimate 24gi-p^ ~ 1. 

Then due to the first term jj ^ -f^-^ which means the maximum of ffq^ in ([12) is less than one. A successful 
inflation in the slow-roll regime restricts the value of the field velocities which may restrict more the above estimation. 
To discuss on the second term let us divide 51(93, x) to its amplitude and functionality as 31(93, x) =| 51 | x/(¥',x) 
such that I (? I is the amplitude of the 51(93, x) and /((p, x) represents its functional form. So the second term can be 
estimated as (by assuming 93 ^ x) 

dip dx 



H \ ' ^^ ' IVP I 2H 



where | gi \ -^ < 1 has to be satisfied. On the other hand one of the Friedmann equations (in the absence of 
the potential) is H^ = ^<p^ + -^X^ ~ IvpLfP' for the zeroth order of 51 -pr- Now if 5^ or ^ have the significant 
amplitude with respect to M then a large amplitude of non-Gaussianity would be expected. This can be realized by 
assuming sharp features in the functionality of 31(93, x) maybe due to a phase transition. 

D. Some Clarifications on Differences with Senatore and Zaldarriaga Q] 

The significant difference is the existence of the terms containing the adiabatic and entropy perturbations themselves 
not just their derivatives. The reason for this difference is in how the effective field theory is constructed in Q. As 
mentioned before in their model the adiabatic mode is borrowed from 2] which satisfies a shift symmetry. Then 
the entropy modes are added and satisfy the shift symmetry too. Consequently, in their formalism they have just 
derivative of perturbations. But in contrast we do not start with distinguishable fields then we do not have any 
difference between the perturbations initially. So by this starting point we had to define the adiabatic and entropy 
perturbations. This is what has been done in this section in details. Now the question is that is there any special 
transformation for adiabatic and entropy perturbation in our model? Yes, it is locally rotated shift transformation 
i.e. 

(5cr ^ (5o- + (ci cos 6* -I- C2 sin 0) (26) 

5s -^ 5s + (ci sin 9 — C2 cos 9) 

where 9 = arctan(x/93). Note that the rotational angle depends on the background fields time evolution. To achieve 
this result, the starting point is the Lagrangian for two fields i.e. the relations pUj) . ([TT|) and (|12p . By looking at 
these relations it is obvious that 5ip ^ 5ip + ci and 5x ^- 5x + C2 is a symmetry of the model where ci and C2 
are two independent arbitrary constants. So due to ((T3)) one can get the above relation ((26)) as the corresponding 
transformation of 5a and 5s. According to this relation the invariant terms corresponding to 5ip and 5x are not 5a 
and 5s but 

5a - e5s -^ 5a - 95s (27) 

5s + 95a -^ 5s + 95a 

which are T.5 and N.5 respectively. Not surprisingly, these terms construct whole Lagrangian in adiabatic and entropy 
perturbations language as seen previously. 

So it seems initially supposed adiabatic perturbation in [J| results in lack of all possible terms in the effective 
Lagrangian. Our proposition to solve this problem is based on the discussion in this subsection. The main building 
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blocks for an effective field theory of multi-field inflation are not Sa and 5s but they are T.6 and N.S. So the most 
general Lagrangian for the perturbations in the multi- field context should be written as[31| 

A/: oc ^ c™„ (f.i\ (N.i\ (28) 

m,n \ / \ / 

for arbitrary Cmn- The above Lagrangian can be considered as the effective field theory for the two-field inflation in 
the language of jj] but with additional terms. Note that the above result can be generalized to multi-field inflation as 

A/:«^c„„,„,,...,„„ffij (Ni.5j (n2-SJ ...(Nn-SJ (29) 

where T and iV^'s are a set of orthonormal vectors for an (N + l)-field model. 

V. CONCLUSIONS 

In this work the effective field theory of multi-field inflation has been studied as a generalization of Weinberg's idea 
[l| for a single field. In this approach the most general Lagrangian is built by using all the covariant terms of the fields. 
Though effectively the terms with higher order derivatives are interested in the higher energy scales. In this work we 
restricted the model to the first correction terms. They results in up to fourth order terms in perturbations. Then due 
to the physical interests we switched to the adiabatic and entropy formalism. It has been shown that generally these 
modes can have different speeds of sound. By considering the non-linear terms we studied the non-Gaussinity in this 
model. It has been shown that the amplitude of non-Gaussianity can be significant when the curvature of the classical 
path in the phase-space becomes large. For example a sharp turn in the classical path can realize it. However it 
seems that existence of the higher order derivative terms in the Lagrangian cannot produce large non-Gaussinity. The 
bottom line for this fact is the strong constraint on the coefficients to keep the effective field theory valid. But there 
is an idea that it is possible to take the higher order correction terms under control automatically, e.g. by Vainshtein 
mechanism. This relaxes the constraint on the coefficient of the correction terms and results in large non-Gaussianity. 
On the other hand the structure of the interacting terms in the Lagrangian predicts the existence of all the shapes 
of non-Gaussinity with the different amplitude for different cases. But the characteristic feature of the model is that 
the non-Gaussinities are correlated. That means if there is a local non-Gaussinity due to the entropy mode then 
certainly there is a non-Gaussinity in adiabatic mode which is equilateral. The amplitude of these different types of 
non-Gaussinities are not independent. 

In contrast to [J|, the adiabatic and entropy perturbations are not distinguishable initially. This fact results in the 
existence of the perturbations as well as their derivatives. In other words the adiabatic and entropy perturbations are 
not invariant under the shift symmetry of original fields. However a combination of them is invariant under such a 

symmetry. These combinations are "5cr — 96s" and "5s -I- 96(j" or in other form T.6 and N.6 respectively. This result 
is important for constructing the effective field theory for multi-field inflation and causes the additional terms with 
respect to what is considered in ^. 
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Appendix A: The Most General Lagrangian for Multi-Field Inflation 

In this appendix we will study the most general form of a Lagrangian including up to the fourth order space-time 
derivatives for multi-field models. It is a straightforward generalization of what has been done in [16] for a single field. 
To continue let us review the blocks 
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which satisiy general covariance. Obviously all the indices can be raised by the metric inverse g^"^. To construct the 
Lagrangian we need to make scalars from the above blocks. The possible terms are listed in |16l] and easily can be 
generalized for more than one field. The most general form of the Lagrangian with the above blocks for a multi field 
model reads as C/ ^ = 

+ c^((/?)(V^i?)(V^(^7) + C5{0)UR + 4((/?)i?^,V^V>/ + ai{0)R^,p„R''''P'' + a2{0)Rf,,R^"' + a^{0)R^ + ai{0)R 

where summation on repeated indices is assumed and = (^i, (^2, •••, '^n for N fields. Note that all of the terms in 
(jAip are not independent up to total derivative terms. In the following we show this fact in details. 

ci(^)(V'^i?)(V^V7,) = V^'{c{{if)RV^,ipI)~c{{0)RU<pI-R{V^'c{{0))V^VI 

= V^{ci{ip)RV^ipi) ~ cii0)Rn^i - R{dci{ip)/d^j)V''ipjV^ipi 

C5i0)nR = V^(c5(</?)V^i?) - (V^C5((^))(V^i?) 

- V^f(c5(^)V^i?) - ((V^C5(^))i?)) + {W^V^'c.mR 



V^( (c5(^)V^i?) - {{VpC5mR)j+Vj{dc5i0)/dipi)V^^i\R 

V^('(c5(^)V^i?) - {{\/^c5{0))R)] + (I^^^M^'/V'^V,/ + (ac5(v?)/5<^/)V^V^^/ ]R 



= b'/{0){Wifi)V^Difj + bl^{0){V^vi)R^,,V''vj + 

- V^[bi-\0){V^^i)U^.!^~{VX-'mi^^Vi)^VJ-hi\0){U^i){Uv.j) 
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difjdipK ) 






Due to these relations, the most general Lagrangian (IA1|) reduces to a Lagrangian with redefined coefficients and 
vanishing fee, 67, &8, &9, &10, &11, C4, C5 and ce- 

In this step let us deduce the equations of motion for the leading terms. To do this let us assume 04 ((^) = — ^, 

^^A kv) — —^S^'^ and b5{(f) — ~MpU{0) by redefinition of the fields without any loss of generality. Now the 
equations of motion for the ipi and y^j, are 

M^U<pi^MldU{0)/d^i, R^, = -{MyMl)S''^V^^iV,ipj~Uiip)g^,. (A2) 

Now by plugging the above relations to the remaining terms of the Lagrangian (jAip . it is easy to show that 



a3(v?)i?' = a3i0)[ -Rf,,p„R^"'P''+4R^,.R^'n ^aii0)R^,p„R^'''P'' + a2i0)R^,R^''. 
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Due to the above relations and again a redefinition of 0, U{(f), 63 "'^^((/3), ai((p) and a2{0) the most general form of 
the Lagrangian becomes 

which exactly reduces to what is deduced in jljj for a single field model. 

Appendix B: Comparison With Gordon et al. [10] 

It is obvious that the LHS of (ITT]) is exactly compatible with equations (47) and (48) in [10|]. The RHS of equation 
(47) in [lO| can be written as 265s + 2(0 — ^9)5s. So to show that our calculation is compatible with this result one 
has to show that 6 + {'iH + 2^)6 + V^s = 0. To show the RHS is also the same let us start by pointing out that 

^=^(^X-X^) = 4(x^^-^^x) -^ 9^-2^{xVv-0V^) + ^{xV^-0V^)^V^s (Bl) 

-^ e + 2-e + K. - ^{xv^ - Wx) = 'iH^iv^v^ - xv^) = -sm qed, 

a a (T 

where in above calculation we used Klein-Gordon equations + 3Hip+V^ = 0, x + 3Hx+V^ = and (j + 3H&+Va = 
due to a^ = LfP'+x^ . To compare the second equation in (jl7|) and (48) in [l^ it should be checked that if the coefficients 
of 5a and 5a are same or not. The former is trivial and the latter results in checking 29^ — —{9 + 3H9 + Va-s)- This 
equality has been already checked in the above. So the result is exactly compatible with |10l |. 
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This term is -l{Va,ySa'^ + V^sSaSs + V^sSs'^). 

Here we do not expand T.5 and N.S since they contain no common terms to factorize. So their expansion may cause just 

messy stuffs without any physical interests. 

Note that we employ the standard definition of Cs . It means we skip the interaction terms between 5a and Ss which exist 

even in quadratic level. 

We here just consider the terms containing time derivatives and not any spatial derivatives. However the procedure is 

same. 

Note that the coefficient of {T.oj in H18|l is 6. But by comparison to H19|) . 2 of 6 appear in definition of the speed of sound, 

Cs, (exactly same as the second line in Q) and what remains is 4(T.o) . 

Note that due to the first equation of motion in (|17p the Ss sources 5a. 

Note that except here in the rest of the paper we assumed that Qi's are constant as a matter of simplification. But here 

we would like to show how the functionality of gt's may affect the final result. 

Note that the most general form of the potential term can be supposed. But as mentioned in '8| , in the slow-roll regime 

there is no interesting non-Gaussianity prediction for single field models. However for multi-field models the potential term 

can result in large non-Gaussianity which considered in 9]. Here, we restrict our calculations to kinetic terms. 

Note that here we just consider the time derivative since in the discussion of this section there is no difference between 

our model and [J] for the terms containing spatial derivatives. This is because the background is not spatial dependent. 

Remember that the angle of rotation just depends on time. 



